Theory of Optical Transmission Through a Near-Field Probe With a Dissipative 

Matter in its Core 
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We develop a theory of light transmission through an aperture-type near-field optical probe with a 
dissipative matter in its semiconducting core described by a complex frequency- dependent dielectric 
function. We evaluate the near-field transmission coefficient of a metallized silicon probe with a 
large taper angle of in the visible and near-infrared wavelength range. It is shown that in this 
spectral range the use of a short silicon probe instead of a glass one allows to achieve a strong (up 
to I0 2 — I0 3 ) enhancement in the transmission efficiency. 



The development of scanning near-field optical mi- 
croscopy (SNOM) has already led to enormous progress 
in studies of different nanoscale phenomena. Among 
other things, this technique is widely used for image for- 
mation of various nanoobjects, single molecule fluores- 
cence detection, and for laser-induced ablation of a sam- 
ple close to the tip apex (see In the recent years 
SNOM methods have been employed in studies of quan- 
tum dots and single- walled carbon nanotubes . The 
near-field optical microscopy has opened a possibility for 
direct observation of fine features in the self-focusing ef- 
fect , for imaging the light propagation in photonic 
crystal waveguides |fj and the electromagnetic local den- 
sity of states of optical corrals Q. 

To improve efficiency of the aperture probes for the 
SNOM technique one needs simultaneously to increase 
their transmittance and spatial resolution capability. In 
the present letter we study theoretically one possible way 
to enhance the optical transmittance through a met- 
allized near-field probe with a subwavelength aperture. 
We suggest to work in the visible region using a short 
probe with a large taper angle with a core consisting of 
a semiconducting matter with a high refractive index n. 
The increase in transmittance with n occurs due to de- 
crease of the light wavelength A c = l-KcjriLO inside the 
core. Consequently, at high n the cutoff effect, which 
strongly reduces the light transmission through the probe 
in the overdamped regime, affects its transmission effi- 
ciency much less than for a glass. However, this gain 
may be counterbalanced by the growth of light absorp- 
tion in the dissipative medium of the core described by 
the imaginary part of the frequency dependent dielectric 
function e(u>). To deal with these two competing tenden- 
cies one needs a detailed theory. 

The purpose of this letter is to develop an analytic 
theory for transmission of visible light through a probe 
with a core made of a dissipative semiconducting mat- 
ter. Here we are encouraged by the transmission SNOM 
experiments with a silicon probe at A = 1.06 fim j3>0|, 
which indicate that such probes are very promising in the 
near-infrared (IR) region. The theoretical support for 



this work comes from a comparative numerical analysis 
of the transmission efficiencies of glass and silicon probes 
at A = 1.3 /im We have to note however, that the 
above theory is restricted by the use of a two-dimensional 
model with a loss-free dielectric core and small taper an- 
gle of 15°. The further experimental success was achieved 
in , where the authors have employed a pyramidal Si 
probe that was entirely coated with a thin metal film to 
increase the transmission efficiency in the near-IR region 
(A = 830 nm). An extremely high throughput (2.3 %) 
was obtained in this experiment with a resolution capa- 
bility about 85 nm. In the IR region the light absorption 
in Si is sufficiently small, but in the visible region the 
imaginary part of its dielectric function rapidly increases 
with a decrease of A. So, dissipation of the electromag- 
netic energy inside the silicon core and frequency dis- 
persion of its dielectric function become important and 
should be taken into account to get an adequate physical 
pattern of light transmission through the probe. 

The three-dimensional theory developed in the present 
letter is valid for all taper angles including large ones 
which are especially suitable for large transmittance 
|12l [l3j . The theory is based on the exact analytic de- 
scription of the conical waveguide eigenmodes inside a 
probe with a dissipative matter in its core and perfectly 
conducting metallic walls. For a loss-free dielectric core, 
similar approach has been developed in our recent work 
|l3l | . Here we consider the case when the semiconducting 
core has a complex dielectric function which is a nec- 
essary feature of the visible region. This generalization 
requires more than a trivial analytic continuation of our 
previous results due to the necessity of describing the 
effects associated with frequency dispersion and light ab- 
sorption inside the core matter. Such effects lead to a 
number of new features in light transmission through 
semiconducting probes, which discuss below. 

We consider here time-harmonic fields inside a cone 
whose core consists of a dissipative medium and whose 
walls are perfectly conducting. In spherical coordinates 
the Helmholtz equation for the Hertz function U of the 
electromagnetic field inside a cone is 
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Here r is the distance from the cone vertex, 9 and (p 
are the polar and azimuthal angles. For a dissipative 
matter, the wave number k is complex: k — wyf/i/c, 
y/ejl = n + in, where c is the speed of light, n and k are 
the refractive index and the attenuation coefficient. As- 
suming the permeability /j, = 1, the real s' = 5Re{e} and 
imaginary e" = 3m {e} parts of the frequency-dependent 
dielectric function e (w) = e' (ui) + is" (uj) are given by 
e' = n 2 - k 2 and e" = 2nn (see Ref.0). 

The relevant solution of Eq. corresponding to the 
standing wave with vanishing amplitude at the cone ver- 
tex (r = 0), has the form 



U 



= K (r) P™ (cos 9) e lmv , K = Crj v (kr) , (2) 

at which the radial dependence 1Z (r) of the Hertz func- 
tion @ is expressed through the spherical Bessel function 
of the first kind j v (z) of a complex argument with the 
index v not equal to an integer. Here C is a constant. 
The dependence on the polar angle 9 is determined by the 
associated Legendre function of the first kind P™ (cos 9) 
with power v and order m (m is an integer). 

For the transverse magnetic (TM) field modes the 
boundary condition at an interface between a core 
medium and a perfectly conducting metallic coating of 
a conical waveguide can be written as P™ (cos6>o) = 0, 
where 9 is the cone half-angle. Each choice of numbers m 
and n in this equation (n denotes the number of its root) 
determines a possible TM mn mode. The eigenvalues v mn 
strongly depend upon the value of 9q such that v mn de- 
creases with an increase of #o- hi the most interesting 
case of the lowest-order TMoi mode, the projections of 
electric E and magnetic H fields onto the corresponding 
axes of spherical coordinates (r, 9, ip), take the form 



E r = -^—}- K(r)P v {cos 9) , E 



dK(r) dP u (cosf 



rdr 



09 



(3) 



1 dW) \Hj\ 



(0) 



167T duj 

To determine the near-field transmission coefficient of a 
truncated conical waveguide we introduce the quantities 



Wp(r) = 2nr 2 / w p (r, 9) sin 



(7) 



which represent the integrals of w r , wg, or w v taken over 
a part of spherical surface lying inside the cone (0 < 
8 < < <p < 2ir) at a given distance r from the 
cone vertex. With the help of Eqs. ©, and I0H0I, 
these integrals can be evaluated explicitly. The resulting 
expressions for W r , Wg, and W v , take the form 



\C\ 2 d{ioe') 



v{v + l)I v j v (n + in) — (8) 
duj L c 



LOT 



(n + in) 



Lur 



Jv+l 



{v + 1 ) j v (n + in) — 
c 



(n + in) 



(9) 



Lor 

{n + in) — .(10) 

c 



Here the angular integral T v is given by 
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The integral energy density W to t — W r + Wg + W v 
can now be evaluated as the sum of Eqs. (|1(J[1 . At 
small distances from the cone vertex (r <^ A c ) it exhibits 
a rapid power fall, Wtot °c (|fc|r) 2,y , with a decrease of 
r. At distances r from the cone vertex much greater 
than the wavelength in the core medium (r 3> A c ), the 
asymptotic expression for Wtot takes the form 
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(4) 



while the H r , Hg, and E v components are equal to zero. 

In a dissipative media with frequency-dependent di- 
electric function e — e'+e" and permeability /i = u'+ifi" , 
the general expressions (see [lj|) for the time- averaged 



densities of the electric, 



w e i 



wg, and magnetic, 



fields are given by 
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Here r K = c/2klu is the attenuation length. It is ev- 
ident that the main feature in the radial dependence 
()12J) is determined by a factor cosh (r/r K ), which re- 
flects the influence of light absorption inside a dissipa- 
tive core of a near-field probe. Another important point 
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is that the integral energy density Wtot(r) exhibits an 
oscillatory behavior at distances far from the cone ver- 
tex r 3> c/u>\n + ik\. These additional oscillations are 
the result of the frequency-dependent dielectric function; 
they are absent if the core is made of a loss-free medium 
(this is the case of a probe with a glass core, e — const 
and k — 0, see 13] ). However, for a lossy matter the am- 
plitudes of the electric (W et cx d(uje') /du, Eqs. ©, ©) 
and magnetic (W m cx |e|, Eq. HlOjl) components of the 
integral energy densities are not equal and can be signif- 
icantly differed from each other. Therefore, the oscilla- 
tions of the electric, W e i, and magnetic, W m , energies do 
not compensate each other far from the cone vertex. Note 
also that in the presence of light absorption there is a sig- 
nificant difference in magnitudes of the energy fluxes as- 
sociated with the incident and the reflected waves. Their 
ratio \S r /Si n \ is equal to exp (— 2r/r K ). So, this difference 
becomes particularly important in the range of r r K . 

Further we evaluate the optical transmittance of a con- 
ical waveguide with a dissipative matter in its core. For 
a probe tapered to a subwavelcngth diameter it is nec- 
essary to distinguish a near-field transmission coefficient 
of a waveguide itself and the result ing transmission co- 
efficient to the far-field zone (see The near-field 
transmission coefficient, T, can be expressed in terms of 
the time-averaged energy densities associated with the 
output and the input fields of the waveguide. For spher- 
ical waves inside a cone, this coefficient can be defined 
as the ratio T = W t °"* /W™ t of the time-averaged energy 
density W t °"* = W to t (zo) at the exit plane z = z of a 
probe integrated over the aperture cross section 2npdp 
with radius a (a — zotan#o) t° the corresponding inte- 
gral energy density W^ t at the waveguide entrance with 
radial coordinate ri n . In a dissipative medium, the lat- 
ter is given by W t l ™ t = aW to t (r.n) (see Q). The fac- 
tor a = [1 + exp (— 2ri„/r K )p 1 shows a fraction of the 
integral energy density at r = ri n , associated with the 
incident wave alone. So, the contribution of the reflected 
wave turns out to be completely removed. 

For the subwavelength aperture 2a <C A c , the basic 
expressions (fHf) — ffTTTf) for W r , We, and W v at the exit 
of a conical waveguide, can be expanded in power se- 
ries of |fc|r out , where r out = a/sin#o is the correspond- 
ing radial coordinate. For the value of Wtot (fin) at the 
waveguide entrance (?-£„ ^> A c ) we use the expression 
W t ot(r) cx cosh (r/r K ), averaged over the fast oscilla- 
tions of Eq. (|12[1 . Then, the resulting expression for the 
near-field transmission coefficient takes the form 

„ / LJ 71 + IK \a \ ,i Tin , ^ 

T cx ' . ' cosh^ 1 — ^- . (13) 

The eigenvalues v = vq\ of the TMoi mode in Eq. i|13|) 
exhibit rapid fall with an increase of the taper angle 
Ooi = 4.083, 2.548, 1.777, and 1 at 6 = tt/6, tt/4, 



7r/3, and n/2, respectively). Thus, Eq. (jT3"j) describes 
well all major features of light transmission through the 
subwavelength aperture in a conical waveguide with a 
dissipative matter in its core. It is seen that the values 
of T are strongly dependent on the ratio a/ A, the taper 
angle 29q, and the refractive index n. Moreover, accord- 
ing to Q13[) the transmission coefficient T is proportional 
to cosh -1 £, where £ = l/r K is the ratio of the length of 
the probe edge / to the attenuation length r K = c/2klo 
(at small a we have ri„ w I), It is clear that the high 
transmission efficiency of a semiconducting probe can be 
achieved in the wavelength region far from the peak in its 
absorption band (k -C n). If additionally r K 3> I, than 
one can put r K — * oo. Then, Eq. I|13|) is reduced to espe- 
cially simple form T cx (w na/c sin 8 o) 2 "*- 00 -*. This is the 
case of a loss-free dielectric core. In the opposite case of 
large losses (r K <C I), the transmission coefficient behaves 
like T cx {una/ csviOq) 2 "^ ^ exp (— £). This reflects the 
strong influence of light absorption on the value of T. 

Now we apply our theory for studies of the trans- 
mittance of the visible and near-IR radiation through 
the aperture-type metallized silicon probe. In the wave- 
length region from 830 nm (fiui = 1.5 ev) down to 400 
nm (fiu> = 3.1 ev), the refractive index n of Si increases 
monotonically from 3.67 to 5.57 and the attenuation co- 
efficient k grows from 0.005 to 0.387. This results in quite 
different influence of light absorption in Si in the near-IR 
and the short-wavelength part of the visible spectrum. 

In Fig. 1 we present the wavelength dependences of 
the near-field transmission coefficient of the metallized 
silicon probe for the most interesting case of large taper 
angle 20q = 90° and for various values of the aperture 
diameter 2a. To demonstrate a dependence of light ab- 
sorption inside the Si core on the length of the probe edge, 
we calculated the transmission coefficient T for various 
values of I. As expected, T is strongly dependent on the 
aperture diameter in full agreement with simple formula 
(|13fl derived in this work. However, the wavelength de- 
pendence, obtained in the present work for the Si probe, 
differs dramatically from the case of a loss-free dielectric 
core (a glass fiber), for which T cx (a/A) 2 "). As is evi- 
dent from Fig. 1, the transmittance of the silicon probe 
strongly varies over the spectrum. If the length of the 
probe edge is not too large (I < 10 /Ltm), the transmis- 
sion coefficient increases first as the wavelength decreases 
from the IR region, reaches its maximum at a definite 
wavelength A max , and then strongly falls at A <^ A max in 
the short-wavelength part of the visible spectrum. The 
position of the maximum A max and the maximal value of 
T max depends on the specific geometrical parameters of 
the probe. 

It is important to stress that this maximum in the 
transmission efficiency of a silicon probe lies at A ~ 550 — 
800 nm (see Fig. 1). This occurs despite the fact that the 
attenuation length r K , associated with the imaginary part 
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FIG. 1: The near-field transmission coefficient, T, of the con- 
ical waveguide (2#o = 90°) with the Si core vs A = 2ttc/ll>. 
Curves 1, 2, 3, and 4 correspond to the aperture diameter 
2a = 100 nm, 70 nm, 50 nm, and 25 nm, respectively. The 
length of the probe edge is I = 2 /an (a), 4 )im (b), 8 /jm (c). 



of the dielectric function of Si, considerably decreases 
in the visible region compared to the near-infrared one. 
For example, at A = 633 nm, 532 nm, and 488 nm, the 
respective values of the attenuation length r K turn out to 
be equal to 2.66 /im, 0.84 /im and 0.45 /im, in contrast 
with 13.53 /*m at A = 830 nm. 

In summary, it follows from our calculations that at 
large taper angles high values of the near-field transmis- 
sion coefficient can be achieved for a passage of visible 
light through Si core of an optical probe. To illustrate the 
enhancement in the trasmittance of Si probes in compar- 
ison with conventional fiber ones we compare the present 
estimates for Si with those obtained for a core with small 
n (glass or S1O2). Although the taper angles of fiber 
probes do not usually exceed 40° which additionally re- 
stricts their efficiency, we use the value 26q = 90° and 



n = 1.55 to make comparison with our recent results 
[l3j . For the probe with the length I — 2 /im and the 
aperture diameter 2a — 50 nm we get Tsi/T g i ass = 2.2, 
14, 45, and 71 for A = 488, 532, 633, and 830 nm, respec- 
tively. For the same parameters, but 29o — 60° we have 
T S i/T g i ass = 43, 240, 800, and 960. According to our the- 
ory, the enhancement occurs as a result of competition 
between two factors: the rise of n and the decrease of the 
attenuation length r K Ijl3|l . As follows from our results, 
the former factor is more important in the most part of 
the visible spectrum, provided the near-field probe length 
is sufficiently short (no more than several /im) , such that 
effects associated with the light absorption are not too 
strong. We also would like to point out that in case of 
an entirely coated probe the surface plasmon-polariton 
propagation along the metallic cladding is likely to fur- 
ther increase the resulting transmittance in accordance 
with the mechanism discussed in Ref. [l2|. 
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